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Abatract 


The  loaa  of  momentum  and  atiffneaa  due  to  ab¬ 
lation  may  significantly  influence  the  vibrations 
of  a  solid  propellant  grain.  This  paper  presents 
an  analytical  study  of  the  axial  shear  vibrations 
of  a  long  hollow  cylinder  that  is  subjected  to 
time  dependent  body  forces  in  the  axial  direction. 
The  outer  surface  of  the  cylinder  is  bonded  to  a 
rigid  case,  and  the  inner  radius  Increases  tnono- 
tonically  with  time.  An  expression  is  determined 
for  the  shear  stress  at  the  bond-interface.  It  is 
shown  that  the  frequency  of  the  shear-bond  stress 
increases,  and  that  its  amplitude  decreases 
towards  burnout  time.  The  shear  stress  is  studied 
for  various  ablation  rates.  Conventional  methods 
of  analysis,  such  as  separation  of  variables  and 
Four ler-Bessel  analysis,  are  not  directly  appli¬ 
cable  in  this  problem,  since  the  boundary  condi¬ 
tions  are  prescribed  on  a  time  dependent  surface. 

A  modified  Fourler-Bessel  mode  is  defined  that 
satlafles  the  boundary  conditions.  By  substitut¬ 
ing  this  mode  into  the  equation  of  motion,  a  solu¬ 
tion  is  obtained  by  asymptotic  methods  in  the 
vicinity  of  the  bond- Interface.  The  analysis  is 
extended  to  Include  the  axial  shear  vlbrationa  of 
an  ablating  viscoelaatic  cylinder.  Viscoelastic¬ 
ity  is  Introduced  by  means  of  the  relaxation 
function  in  shear. 
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I.  Introduction 

The  designing  of  solid  propellant  rocket 
motors  requires  a  consideration  of  the  dynamic 
response  of  a  propellant-casing  system  to  time- 
dependent  body  forces.  The  body  forces  may  be 
associated  with  spinning  motions,  axial  sccelera- 
tions,  etc.  This  paper  focusea  upon  the  axially- 
sysnetrlc  dynamic  response  due  to  axial  accelera¬ 
tion,  such  as  occurs  in  accelerated  flight. 

Special  attention  is  devoted  to  the  shear  bond 
stress  at  the  propellant-casing  Interface. 

The  free  and  forced  vibrations  of  encased 
elastic  cylinders  have  been  considered  for  various 
types  of  surface  constraints.  The  dynamics 

of  encased  viscoelastic  cylinders  have  also  bean 
studied.  These  paperu  present  results  for 
various  values  of  the  ratio  of  the  Internal  and 
external  radii  of  the  propellant  cylinder.  In 
this  way  one  obtains  quasl-statlc  information  on 
the  Influence  of  an  increasing  inner  radius.  It 
was  only  recently  that  the  Influence  of  continuous 
ablation  on  the  axially- symmetric  plane  strain 
vibration  was  studied  for  both  an  encased  elaetlc* 
and  an  encased  viscoelastic  cylinder.?  In  the 
present  paper  an  attempt  is  made  to  analyxe  the 
Influence  of  continuous  ablation  on  the  axial 
response  of  a  burning  rocket.  For  siaq>llcity  we 
consider  a  case-bonded  grain  of  infinite  length. 

To  simplify  the  analysis  the  presence  of  any 
star  points  is  neglected,  and  the  propellant  is 
represented  as  a  thick-walled  cylinder.  Any  star 
point  material  nay  be  accounted  for,  however,  by 
taking  the  inner  radius  of  the  propellant  cylinder 
as  the  radius  the  cylinder  would  have  if  the  star 
point  material  were  uniformly  distributed  around 
the  inner  surface.  It  is  assumed  that  the  stiff¬ 
ness  of  the  propellant  is  so  small  compared  to  the 
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stiffness  of  the  casing,  that  the  latter  may  be 
aaaumed  as  Infinite.  In  a  study  of  a  related 
free-vlbratlon  problem  by  Baltrukonie^  It  was 
shown  that  for  axially  symmetric  shear  deforma¬ 
tions  thle  assumption  la  acceptable  for  the  stiff¬ 
ness  ratios  that  era  commonly  encountered  In  solid 
propellant  motors. 

The  paper  la  prlamrlly  concerned  with  the  In¬ 
fluence  of  grain  ablation  on  the  shear  bond 
•trees.  It  la  shown  that  the  loss  of  mass  and 
stiffness  of  a  burning  grain  significantly  affects 
the  magnitude  and  the  frequency  of  the  shear  bond 
stress,  fines  most  propellants  are  viscoelastic 
In  shear  the  Influence  of  viscoelastic  damping  la 
Included  In  the  third  part  of  the  paper. 

The  dynamic  solutions  that  are  presented  In 
this  paper  are  discussed  In  relation  to  the  solu¬ 
tions  of  the  analogous  quaal-atatlc  problems.  The 
quasl-statlc  solutions  for  the  elastic  cylinder 
are  simple.  A  quaal-atatlc  solution  for  the 
ablating  viscoelastic  cylinder  Is  presented  by 
Lindsey  and  Williams. * 

The  paper  consists  of  three  parts  which  are 

concerned  with  the  axial  shear  vibrations  of  (1) 

an  encased  elastic  cylinder  of  constant  Inner 

radius  a  ,  (2)  an  encased  elastic  cylinder  of 
o 

monotonlcally  Increasing  radius  a(t)  ,  and  (3)  an 
encased  viscoelastic  cylinder  of  monotonlcally 
Increasing  radius  a(t)  .  Conventional  methods  of 
analysis,  such  as  separation  of  variables  and 
Fourier-Besael  analysis,  are  not  directly  appli¬ 
cable  In  the  last  two  problems,  since  the  boundary 
conditions  are  prescribed  on  a  time  dependent  sur¬ 
face.  A  modified  Fourler-Beeael  mode  la  defined 
that  satisfies  the  boundary  conditions.  By  sub¬ 
stituting  this  mode  Into  the  equation  of  motion,  a 
solution  Is  obtained  by  asjrmptotlc  methods  In  the 
vicinity  of  the  bond- Inter  face.  The  method  thus 
allows  ue  to  determine  the  shear-bond  stress. 


u  •  u  *0 
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The  cylinder  is  Infinitely  long  and  the  system  Is 
antl-synmetrlc  relative  to  any  plane  perpendicular 
to  the  i-axls.  The  stress  component  o  ,  which  la 

a  sysmsetric  quantity,  must  then  vanish.  The 
resMinlng  stress  and  displacement  cosmonauts  are 
Independent  of  c,  and  equations  (2)  reduce  to  two 
uncoupled  equations  of  motion.  The  equation  of 
motion  In  radial  direction,  subject  to  the  bound¬ 
ary  conditions  (3)  and  quiescent  Initial  con¬ 
ditions,  yields  only  the  trivial  solutions 
°rr  "  cee  "  0  •  The  only  remaining  equation  con¬ 
stitutes  the  governing  equation  of  the  present 
problem. 

r  h  <r°r«)  4  p,i(t)  "  p“t  (4) 

The  solution  of  Eq.  (4)  la  subject  to  the  boundary 
conditions 

at  r  -  a(t)  org  ■  0  (5a) 

at  r  ■  b  u  •  0  (5b) 

m 

Assuming  an  initially  undisturbed  cylinder  the 
Initial  conditions  are  expressed  as 

u  (r,t)  -  u  (r,t)  -  0  for  t  <  0  (6) 

s  s 

We  shall  determine  solutions  for  an  elastic 
cylinder  as  well  as  for  a  viscoelastic  cylinder. 

By  eliminating  the  stress  from  Eq.  (4),  the 
equation  for  the  axial  displacement  ug(r,t)  of  the 

elastic  cylinder  is  obtained  aa 

7  37  G  Sr)  +  «.<*>  ■  <7> 


11.  The  Equations  of  Motion 


In  a  solid  body  the  stress  tensor  satis¬ 
fice  the  equation  of  motion 


“ij.j  *  wi(V°  '  p”i 


where  F£(Xj,t)  is  a  body  force  per  unit  mass.  A 

uniform  body-force  distribution  Is  assumed,  with  a 
component  In  s-dlrectlon  only.  Since  solutions  of 
Eq.  (1)  are  sought  for  a  circular-cylindrical  body 
subjected  to  axially  aysmetrlc  loading  and  bound¬ 
ary  conditions  Eq.  (1)  reduces  to 


do  So  a 
_ r *  ,  _ at  .  r*  .  n_ 

TT +-3T+T  +  P,s<t) 


(2a) 

(2b) 


The  analysis  la  restricted  to  an  Infinitely  long 
hollow  cylinder  that  la  rigidly  encased  at  the 
outer  surface  and  stress  free  at  the  time  depend¬ 
ent  Inner  surface.  The  boundary  conditions  for 
this  configuration  era 

r  -  a(t)  or  -  org  -  0  (3a) 


If  the  cylinder  Is  linearly  viscoelastic  the  rela¬ 
tion  between  shear  stress  and  shear  strain  is  ex¬ 
pressed  In  the  form 

t 

<*rf  •  J  G(t-s)d(8ug/dr)  (8) 


where  C(t)  is  the  relaxation  function  In  shear. 

The  equation  of  awtlon  for  a  viscoelastic  cylinder 
Is  obtained  by  substitution  of  Eq.  (8)  Into  Eq. 
(4). 

For  convenience  the  following  dimensionless 
quantities  are  Introduced 


ug/b 

T  -  t/t£ 

(9a) 

r/b 

N(T)  •  tfaFg(T)/b 

(9b) 

(C/p)*tf/b 

B  »  a  /b 
o 

(9c) 

In  Eq.  (0)  t  Is  the  total  burning  time.  «e  also 
define  1 

•(t)  -  a  «(t)  ,  where  1  *  or(T)  a  (b/a  >  (9d) 

o  o 

The  governing  equation  for  the  elastic  cylinder, 


* 
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Eq.  (7),  can  now  be  rewritten  as 


(10) 


The  governing  equetion  for  the  viscoelastic  cylin¬ 
der  la 


1  (0  B )Y  (f.  )  -  J  )Y  (.»  ') 

2  onon  n  n 


•f 


2 _ 1__ 

na  U  • 
n 


(20) 


f“  ||  [»  J  l(T-*)d(SW/aEr  ♦  N(T)  ■  “y  (11) 


where  the  following  aubacltutlon  for  the  relax¬ 
ation  function  waa  aade 

C(T)  -  Cg(t)  (12) 

The  Initial  condltlona  on  the  dlaenelonleee  die- 
placement  are 

W(R,0)  -  W(R,0)  -  0  (13) 

In  tenia  of  the  dlmenalonleaa  dleplacement  the 
boundary  condltlona  are  expreeaed  aa 

at  R  -  Per  (t)  du/dt  -  0  (14a) 

at  R  -  1  H  -  0  (14b) 

III.  Elaatlc  Crain  of  Conatant  Inner  Radiua 

For  a  cylinder  of  conatant  Inner  radiua  the 
boundary  condltlona  eiaqilify  to 

at  R  -  B  dW/dR  -  0  (15a) 

at  R  -  1  W  -  0  (15b) 

If  a  auddenly  applied  uniform  body-force  dlstrlbu- 
tlon  la  conaldered  H(t)  ■  NH(t).  The  forced 
vibration  problem  la  defined  by  the  governing 
equation  (10)  and  the  Initial  (13)  and  boundary 
condltlona  (15).  The  aolution  to  thla  problem  la 
aought  In  the  form  of  a  Fourler-Beaael  exp  ana  Ion 
in  terma  of  the  correapondlng  free  vibration 
modea.  The  free  vibration  problem  correapondlng 
to  Eqs.  (10)  and  (15)  waa  conaldered  by  Baltru- 
konia.  The  modea  of  free  vibration  are 


The  aolution  of  the  forced  vibration  problem  la 
written  in  the  form  of  a  Four ler-Benel  expanalon 
aa 

<r 

W(K,t)  -  7  T  (T)U(O  R)  (21) 

L-  n  n 

n»0 

where  the  sum  la  taken  over  the  poaltlve  roota  of 
the  frequency  equation  (17).  In  Eq.  (21)  T  (T) 

are  time  dependent  Fourler-Beaael  coefficients 
which  are,  aa  yet,  unknown  funetlona  of  time.  The 
forcing  function  N(t)  la  also  expanded  In  terms  of 
U(fi  ,R).  For  the  suddenly  applied  uniform  body- 

force  distribution  the  Fourler-Beaael  coefficients 

A  of  NR(t)  are  obtained  as 
n 

An  ■  TTRj/^P)/  [V<nn®>  *  J0,(nn>]  (22) 

By  substituting  the  Fourler-Bessel  expansions  of 
W(R,t)  and  NH(t)  into  the  governing  equation 
(10),  and  by  aasumlng  that  term  by  term  differen¬ 
tiation  la  allowable,  it  is  found  that  T  (t) 

n 

satisfies  the  ordinary  differential  equation 

+  1*nB*T  *  A„  ‘  0  <23> 
non 

The  aolution  of  Rq.  (23)  that  aatlsfles  the 
Initial  conditions  (13)  is 

T  (T)  -  -  (An/nV>cos(0  qr)  4  A/flV  (24) 
n  n  n  n  n  n 

By  substituting  Eq.  (24)  and  Rq.  (17)  into  Eq. 

(21)  the  displacement  solution  of  the  problem  at 
hand  la  obtained  as 

n4  W(1’T>  " 


W  -  U(nn,R)e 


(16) 


where 


U(fl.R)  -  JCQ  R)Y  (0  )  -  J  (0  )T  (O  R)  (17) 
n  onon  onon 

The  elgenf  requeue  las  0^  are  the  poaltlve  roots  of 
the  frequency  equation 


y  J.a(nD6)c..(n„qT)[J<,(n|i«)To(n||) . J„(n.)ro<Q.»)] 

-°  n„,[j.*<n„e>-j„*‘n.>] 

;  Jia(nne)[J0(n,t)t0(n,)  -  vwn.»>!  (JJ) 

«-0  nI,,[J',(nn#)  ■  Jo*<nn>] 


Ji(0  P)Yo(0  )  -  J  (nn>M°  *)  -  0  (18) 

non  o  n  o 

It  is  not  difficult  to  show  that  tht  funetlona 

U(0  ,R)  fora  an  orthogonal  lyittv  on  [8,1 3,  with 
n 

weight  R. 


We  alao  consider  the  static  solution  of  the 
present  problem.  The  static  solution  W*(R)  satis¬ 
fies  the  equation 


d^W*  1  dW*  N 
d?“  R  dR  “  ’  7 


H(T) 


(26) 


1 

f  R  U(Q  ,R)U(0  ,R)dR  -  B  6 
v  n  ■  tm 

8 


(19) 


The  solution  of  Eq.  (26)  that  satisfies  the  bound¬ 
ary  conditions  (15)  Is  easily  obtained  aa 

W*(R)  -  (M/2q*)[p*f n(R)  -  (1/2)(H*-  1)]l(T)  (27) 


where  6  la  the  Kronecker  delta,  and 
ran 


It  can  now  be  shown  that  the  second  expansion  in 


Bq.  (23)  le  the  Fourler-Beseel  nrlii  with  respect 
to  Che  eyetea  (17)  of  (q'/^M)  tines  the  etstic 
solution  V*(t).  In  the  ueuel  feshlon  the  solution 
of  Che  dynaalc  problea,  Iq.  (23),  thus  consists  of 
e  vibration  shout  the  eteclc  equlllbrlua  position. 
Vo  write 

V(l,t)  .  V*(B)  ♦  W(B,t)  (28) 

Of  pertlculer  Interest  In  the  present  problea 
le  the  dynaalc  overstress  et  the  cylinder-casing 
let erf see.  Let  the  dlaenslonlese  sheer  stress  * 
he  defined  ee 


fines  the  boundary  conditions  onW*(l,f),  Bq. 
(14),  ere  prescribed  on  e  tlae  dependent  boundery, 
the  conventional  aethods  of  separation  of  vari¬ 
ables  and  Fourier  analysis  break  down.  In  this 
paper  It  le  proposed  to  seek  a  solution  of  Bq. 

(33)  In  terns  of  erosive  free  vibration  nodes  that 
are  defined  as 

».  [*.<*>•*]  ■ J.  [vt>*3 T.  [vt>1  - 

- J.  [v>]  *.  [vt>»]  <”> 


L  •  a  /pp  b 
rs  ra  a 


(2*> 


The  dynaalc  overetreee  can  then  be  written  as 
Cr,  *  /■>*/»»  (30) 


Bq.  (37)  represents  the  aodes  of  free  vibration 
of  the  constant-radius  problea,  Bq.  (17),  If  the 
tlae  dependent  functions  X  (t)  are  replaced  by 

the  constants  0 i  ,  It  Is  noted  that  Bq.  (37) 

D 

autoaatlcally  satisfies 


By  eaploylag  Bq.  (23)  the  shear-bond  overstrese  at 
I  ■  1  li  date rained  as 


Vn  [Xn(T)»R]  “  0  for  I  -  1  (38) 


Z  Ji*(Cl  ®)coe(0  qT) 

£  •  2  )  - - - 2 - - - rn 

rB  io  o;[v<p.d-j0,<p.>] 

In  Bq.  (31)  the  following  identity  froa 
of  Bessel  functions  was  ueed^O 

Ji(«)T0(s)  -  J0(s)Yt(s)  -  2/ns 

IT.  The  Ablatins  Blaatic  Crain 


(31) 

the  theory 

(32) 


We  now  require 

SVn  [Xn(T)’R]  /  **  *  0  for  *  "  *<T>  <39) 
The  condition  (39)  yields  an  equation  for  X  (t) 

U 

J,  [x„(T)to(T)]  To  [l.(T>]  - 
*  Jo  [Xn(T)]  Tl  [Xn(T)0Qr(T)]  “  0  (40) 


In  the  previous  section  it  is  shown  that  for 
the  non-ablating  cylinder  the  dleplaeaaent  W(B,t) 
can  be  expressed  as  the  sua  of  the  static  solution 
V*(K)  and  a  periodic  function  9(1, t).  It  le  now 
assuaed  that  for  the  ablating  cylinder  the  dls- 
placaaent  W*(1,t)  can  also  be  expressed  la  the 
fora 

wV.t)  -  /‘(B.t)  ♦  9*<1,t)  (33) 

In  Bq.  (33)  ¥**(!, t)  le  the  solution  of  the  quaal- 
statlc  problea,  l.e.  the  solution  of  Bq.  (26)  that 
satisfies  the  boundary  conditions  (14).  If  tha 
body-force  eystea  le  suddenly  applied,  we  easily 
obtain 

«**(!, t)  -  (R/2q* )[8^tf* (v)la(K)  - 

-  (1/2)(lf-  1  )]■(▼)  (34) 

By  substitution  of  Bq.  (33)  for  la  Bq. 

(10)  It  is  found  that  for  t  >  0  the  unknown 
function  9*(1,t)  satisfies  the  equation 


1  89* 


] 


»*g*  t 


05) 


The  Initial  conditions  on  5*(1,t)  at  t  »  0+  are 
obtelned  free  Eqs.  (33)  and  (34)  as 


W*  (1,0+)  -  -  (V/2/)Cf,fn(B)-  (1/2)(I*-1)]  (36a) 

9*  (1,0+)  -  -  (V/2q*)C28*  o-(0+)  fn(l)3  (36b) 

The  function  9*(l,T)  aust  also  satisfy  the  bound¬ 
ary  conditions  (14). 


Equation  (40)  has  the  appearance  of  a  frequency 
equation,  but  the  "frequencies"  are  tlae  dependent 
functions  Xb(t).  On  a  constant  boundary  K  ■  8  the 

frequency  equation  (40)  reduces,  of  course  ,  tr  Bq. 

(18).  The  frequency  equation  (18)  of  the  constant 

inner  radius  problea  was  considered  by  Baltru- 

konli,1  who  tabulated  the  elgenfrequencles 

fi  -  n  (1  -  8  )  for  various  values  of  8.  Tor  the 
n  n  ^ 

first  three  aodes  the  tabulated  H  are  shown  as 

n 

functions  of  6  in  Tig.  1 .  It  can  be  noted  that  ft 

n 

Increases  for  increasing  8,  where  0  approaches 

tx 

Infinity  as  8  approaches  unity.  The  curves  of 
Fig.  1  are  used  to  determine  the  solutions  of  Eq. 
(40).  For  a  certain  tine  t  -  t4  the  solution 

Z  (Tj)  is  obtained  as  [1-8»(t.  tines  the 
n 

ordinate  which  corresponds  ue  abiclssus 
0or(Tl).  In  feet,  since  tht  --eves  in  Fig.  1  are 
aiaost  straight  lines  we  are  justified  in  employ¬ 
ing  a  linear  approxiaatlon.  More  specifically,  we 
can  write  for  8  *  0.333 

MT)  ■  [2.0954-  ,53Bor(T)]  /  [l-8or(T)]  (41a) 

X,(T)  -  [4.9577-  .2488o'(t)]  /  [l-Bar(T)]  (41b) 

Xi(T)  -  [8.0033-  .l5l8or(T)]  /  [l-Bor(T)]  (41c) 

The  expressions  (41)  are  valid  for  8»(t)  <  0.99, 
l.e.  the  present  considerations  are  pertinent  till 
just  before  burnout. 

A  solution  of  Eq.  (35)  is  now  attenpted  in  the 
fora 


0(*/ q)  (48s) 


^(M)  .  /  S  (T)V  (R,T) 
n 


(42) 


S  <t) 
ni' 


r 

X  (  s  )  d  &  ♦ 

^  n 


It  It  noticed  that  the  tenu  of  V^R.t)  are  not 

classical  separation  of  variables  solutions,  sines 

both  S  (t)  and  V  (R,t)  are  functions  of  the  dlmsn- 
n  n 

slonless  time  t.  Nevertheless  Eq.  (42)  Is  eubstl- 
tuted  In  Eq.  (35).  By  assuming  that  the  nth  tens 
of  the  expansion  satisfies  Eq.  (35),  the  substitu¬ 
tion  results  In 

-  9*X*  (t)S  (t)V  (R,T)  - 
n  n  n 

•  !?r  Lv,>V*-T>]  ♦  (43> 


S  (T) 
na'  ' 


T 

X  (s)ds  4  o(3  / q)  (48b) 

o 


It  Is  easy  to  check  that  the  Wrcnsklan  of  S  and 

ni 

S  Is  -21q  .  By  using  a  well  known  theorem^  the 
na 

approximate  solution  of  Eq.  (47)  is  obtained  as 


S„(T) 


sin 


♦ 


The  last  tern  In  Eq.  (43)  Is  yielded  by  the 
Fourler-Bessel  expansion  of  o  w**/frr8  In  terns  of 
the  modes  of  free  vibration,  Eq.  (17),  of  the 
cylinder  of  constant  Inner  radius.  The  function 
Q(t)  and  the  Fourler-Bessel  coefficient  D  are 
derived  as  n 

QCO  -  (N03 /qa )C«(T)8  +o(t)S(T)]  ,  (44) 


♦  PCO 


(49a) 


where 

P(T>  ■  in  I  [*> (T>*» <5)  -  *•  <T>*i  (5>]q(?)«j?  (49b) 

s  o 


and 

V  <"„«>  /  [V  <n0«)-J0*  <«,>]  <«> 

In  an  attempt  to  eliminate  B  from  Eq.  (43)  we 

divide  through  by  V  (B,t).  As  was  expected 
o 

V  (R,t)  does  not  simply  cancel  out,  and  an  exact 
n 

solution  of  the  type  (42)  is  apparently  not  pos¬ 
sible.  He  ahall  now,  however,  restrict  the  an¬ 
alysis  to  the  vicinity  of  R  ■  1  .  Using 
L 'Hospital' a  rule  it  can  be  shown  that 


A  simplification  of  Bq.  (49b)  is  achieved  by 
introducing  the  new  variable 


(30) 


Bq.  (49b)  can  then  be  rewritten  as 

P(T)  -  -  1  [xo(t)]  *1*  {elqv'fS(v)e'tqvdv}  (5ia) 


11m  tt/V  J  -  0 

E-1 

(46a) 

11m  Cv  /V  ]  -  0 

R-1  “  15 

(46b) 

lim  [U(fi  ,R)/V J  -  1 

irl  “  B 

(46c) 

In  view  of  the  above  limits  we  may  write  in  the 
vicinity  of  I  ■  1 ,  that  is,  near  the  bond  inter¬ 
face, 

t(f)  ♦  <«X<T>8n<T>  "  ’  Q<T>Dn  <*7> 

n  n  n  ii 

The  function  B  (t)  is  defined  by  Bq.  (41)  for  the 

D 

first  three  modes.  It  is  apparent  that  in  general 
Eq.  (47)  cannot  be  solved  exactly  for  arbitrary 
ar(T),  But  an  advantageous  feature  of  Bq.  (47)  is 
that  the  real  parameter  q  Is  large.  An  approxi¬ 
mate  method  for  solving  Bq.  (47)  for  large  q  is 
available. 

We  shall  first  consider  the  homogeneous 
equation.  An  asymptotic  solution  of  the  homogene¬ 
ous  equation  can  be  determined  by  Horn's  method.** 
The  solution  which  la  obtained  by  Horn's  method  is 
in  the  form  of  a  eat  lee  of  descending  powers  of  q, 
which  are  asymptotic  to  exact  solutions  of  the 
differential  equation.  By  using  Horn's  method  the 
two  Independent  complementary  solutions  of  Bq. 

(47)  are  obtained  as^ 


In  which 

*1  -  J  *n(«>di 

o 

tad 


(51b) 


♦  (*)  -  Xn  [?(v)]"*/#Q  [?(v>]  (51c) 

In  evaluating  the  integral  in  Bq.  (51a)  we 
again  take  advantage  of  the  fact  that  the  exponen¬ 
tial  exp(-iqv)  contains  the  very  large  real 
parameter  q.  Integrals  containing  such  an  ex¬ 
ponential  are  suited  for  evaluation  by  the  method 
of  stationary  phase. **  The  method  of  stationary 
phase  entails  the  replacement  of  the  real 
Integral  by  a  contour  integral  along  the  lines 
v  ■  -iw  and  v  -  v^-lw.  fines  q  is  a  large  real 
parameter  the  exponential  exp(-qw)  dies  out  very 
rapidly  and  the  main  contributions  to  the 
integral  come  from  near  v  ■  0  and  v  •  vj .  The 
approximate  evaluation  of  Bq.  (51a)  yields 

«*>•  ♦q"*XB(T)"'xn(04)*#,,Q(0*)cos[qJxo(e)ds]  - 

o 

-  q-Xn(T)-Q(T)  4  0(1 /q)  (52) 

The  constants  C*  and  in  Bq.  (49a)  are 
n  o 


T 


V 


deterwlned  fro*  the  initial  conditions  on  t  (t). 

.  0 

It  It  Man  fron  Eq.  (49b)  that  9(0)  •  9(0)  ■  0  . 
Tha  Initial  value  of  5*(B,t),  Eq.  (36),  la  than 

^(t,°)  •  T  [xB(0)]  yi.O)  (33a) 


Also,  by  naglaetlof  terms  of  order  0(1 /q),  we 
derive 


5*«.o)  -  ^  c' 

n 


(1.0) 


(33b) 


Tha  Initial  conditions  on  5*(1,t)  art  praacrlbad 
by  Eq.  (36).  Tha  expressions  In  Eqa.  (36a)  and 
(36b)  ara  In  terns  of  B,  and  thty  ean  be  expanded 
In  Fourler-Baaaal  aarlas  of  U(Ob,B).  It  mi 

obatrvad  bafora  that  0(0.,!)  ■  V  (1,0)  ,  and  tha 

conatanta  C”  and  C*  can  ba  obtained  by  equating 

D  II 


Eqt.  (33a)  and  (33b)  to  tha  Fourier- Baaaal  expan¬ 
el  one  of  Eqa.  (36a)  and  (36b)  respectively.  In 
thle  way  tha  conatant  <r  la  obtained  aa 


37<v>  -  J0'<pm) 


(34) 


Tha  conatant  la  alao  anally  obtained.  It  turne 

out  that  C1  la  of  ordar  0(1  /q)  ae  conger  ad  to  , 
,  B  ■ 
and  C  la  there  fora  naglactad.  lee  paction  of  Eqa. 

(32),  (44)  and  (34)  shove  that  9(t)  la  of  ordar 

0(1  /q* )  aa  coopered  to  C*  ,  thua  9(t)  la  alao 

n 

naglactad.  In  tha  vicinity  of  I  >  1  tha  function 
/(1,t)  nay  than  ba  written  aa 


w\b,t).  -nf  X  C(T)VB(B,T)eoa[q  J  Xn(e)ds]  (55) 
q  n-0  o 

where 


C(t) 


«„<’>'*  Vp.»> 

J'*<na8>  *  Jo'(V 


(58) 


The  dynanlc  ovaratraaa  at  tha  ahaar-bond 
Interface  I  ■  1  li  finally  obtained  aa 

E*(B,t)  -  (q*/»)&3*/&B 

m 

-  2  X  C(t)coe[q  J  Xn(e)de]  (57) 
n-0  o 

Eq.  (57)  properly  raducae  to  Eq.  (31)  for  a  non- 

ablating  cylinder  whan  X  (t)  -  X  (0)  •  0  . 

o  n  d 

V.  The  Ablating  Vlecoalaetlc  Crain 

The  equation  that  governe  the  dynanlc  response 
of  a  vlecoalaetlc  core  la  derived  In  taction  II, 
eee  Eq.  (11).  ly  renovlng  the  dlecontlnulty  at 
t  •  0,  end  by  Integrating  by  parte  Eq.  (11)  le 


rewritten  ee 


si* 

B  51 


.*  Is  1  *  2“ 


d 

Si 


5r  d,J 


daW 


(58) 


The  dlnenelonleee  perenater  q  le  defined  by  Eq. 
(9c),  where  C  now  denotee  the  glaeay  ahear  nodulue 
of  the  vlecoalaetlc  enter lal.  It  le  noted  that 
the  body-force  dletrlbutlon  le  unlforn  and  audden- 
ly  applied. 


A  linear  vlecoalaetlc  notarial  can  In  general 
be  characterised  by  a  discrete  epectrun  of  relax¬ 
ation  tinea  t  .  The  relaxation  function  G(t)  nay 

then  be  axpreeeed  aa 
t 

C(T)  -  CE  ♦  X  Cl  **  <*T/Ti>  (59) 

1-1 


In  Eq.  (59)  G^  denotee  the  rubbery  ahear  nodulue. 

By  canparleon  with  Eq.  (12)  the  function  g(T)  is 
obtained  aa 


t 

g(f)  -  (C^/C)  ♦  X  (Ci/C>  «*P  (60) 

1-1 


In  thle  paper  wa  shall  coaslder  vlecoalaetlc 
natarlale  whoaa  ralaxatlon  functions  In  shear  show 
a  rapid  decrease  for  very  short  tines,  and  then  a 
gradual  decrease  to  the  rubbery  shear  nodulue.  In 
teraa  of  tha  ralaxatlon  spactrun  thle  naans  that 
the  disersta  ralaxatlon  spactrun  consists  of  a 
nunber  of  very  snail  relaxation  tines  (1  <  1  *  l x  ) 
and  a  nunber  of  larger  ralaxatlon  tines 
(Jtj  <  1  s  1).  Mora  specifically  we  oesuna  that 
the  vlecoalaetlc  notarial  con  be  characterised  by 
a  discrete  relaxation  spactrun  such  that 

(1/Tt)»qfi|  for  1  *  1  <  1,  (61a) 

and 

(1/Tj)  «  for  <  1  *  1  (61b) 

In  Eqa.  (61a)  and  (61b)  qO*  Is  the  first  natural 
frequency  of  an  elastic  core  with  the  glassy 
nodulue  ae  shear  nodulue.  The  function  g(T),  Eq. 
(60),  le  now  rewritten  as 

B(t)  -  (G^/C)  ♦  gg(T)  ♦  fc(T)  (62) 

In  Eq.  (62)  gg(T)  la  the  suanatlon  of  exponentials 
over  the  vary  short  relaxation  tinea,  Eq.  (61a), 
and  gs(t)  covers  tha  relaxation  tinea  defined  by 
Eq.  (61b).  The  expression  (62)  ^differentiated 
end  subsequently  substituted  for  g  (t-s)  In  Eq. 
(36).  In  view  of  the^etlpulatlon  (61a)  the 
Integral  containing  gg(T-a)  can  be  elnpllfled. 


I 


I 


1 


fcq.  (38)  can  t lien  be  rewritten  as 


vr  k  [R  S  ] +  r-  k  r.R  J  &,/(T'#)  I i d#] + 


The  vlacoalastlc  material  can  also  be  charac¬ 
terized  by  a  discrete  spectrum  of  retardation 
times.  This  characterization  implies  r  at  the 
creep  function  may  be  written  as 


+  NH  (t)  - 


a2w 

dx3 


(64) 


where 


pa  •  q3  ['  *  Y  <Gi/C>] 

1-1 


(65) 


We  shall  again  seek  a  displacement  solution  in 
the  form  of  an  gsclllatlon  about  che  quasi-static 
displacement  W^  (R,t). 


wJ(R.T)  •  wJ*(E.t)  +  wJ(E,t) 


(66) 


By  substitution  of  Eq.  (66)  for  W(E,x)  in  Iq.  (64) 
the  governing  equation  for  W*(E,t)  is  obtained  as 

if  §i  [R  w\  +  e"  5e  [r  j+fc<T-*>  ST  di]  " 

o 


(67) 


The  initial  condltlona  on  W^(E,T)  follow  from  the 

values  of  the  quasi-static  dlsplacanmnt  and  its 
time  derivative  at  t  •  0*.  The  quasi-static 
viscoelastic  displacement  is  discussed  nest. 

The  stress  for  the  quasi-static  problem  is 

governed  by  Eq.  (4)  if  the  right-hand  side  of  that 
equation  is  zero.  The  dimensionless  quasi-static 
stress  that  also  satisfies  the  boundary  conditions 
(14)  is  easily  obtained  as 


Erz(R*T)  *  <1/2)  [83°'8(T)/R  *  *]  «(T> 


(68) 


It  is  noted  that  Eq.  (68)  does  not  contain  an 
elastic  constant  and  che  stress  solution  is  thus 
independent  of  the  material  behavior.  It  follcars 
that  Eq.  (68)  is  also  the  quasi-static  stress  for 
a  viscoelastic  grain.  It  was  pointed  out  else¬ 
where  that  the  quasi-static  viscoelastic  strain  is 
then  formally  obtained  by  means  of  the  creep 
func  :lon  as 


* 

aw$*/*E  -  j  D(T-s)dE(E,a) 


(68) 


In  Eq.  (69)  D(t)  is  the  dimensionless  creep  func¬ 
tion  .  Integration  with  resoect  to  E  yields  the 
quasi-static  displacement  W”*(g,T). 

T 

j  D(T-a)djyV(s)  fn(E>-  y(«*-1)] 


D(f)  - 


(N/qS  ) 


1  + 


V1 


-x/i 


-  e 


1) 


(71) 


In  Eq.  (71)  d^  are  dimensionless  constants  and 

are  retardation  times.  The  assumption  that  the 
relaxation  spectrum  consists  of  a  number  of  very 
small  relaxation  times  and  a  number  of  larger 
relaxation  times  implies  a  similar  distribution  of 
retardation  times  in  the  retardation  spectrum. 

Let  Dg(T)  cover  the  exponentials  with  retardation 
tlsws  much  larger  than  the  period  associated  with 
the  first  natural  frequency  of  shear  vibrations. 
Following  the  same  procedure  which  yielded  Eq. 

(67) ,  the  quasi-static  viscoelastic  solution  can 
be  rewritten  as 

tfj*  -  [e*«*(T)  in(E)  -  ^i*-1)]H(T)  ♦ 


♦  J  Db(T-a)dE(E,e) 


(72) 


where 


M  •  «  ♦ 


l 


(73) 


The  initial  conditions  on  are  now  determined 

from  Iqs.  (66)  and  (72).  It  is  noted  that  by  this 
procedure  the  influence  of  the  very  short  retarda¬ 
tion  times  is  expressed  in  the  initial  conditions 

on  w!j  rather  than  in  the  forcing  term  on 

the  right-hand  side  of  Eq.  (67).  As  in  the 
elastic  problem  (section  IV)  the  particular  solu¬ 
tion  dueato  the  forcing  function  is  of  order 
0(1/T^aq  )  as  coopered  to  the  solution  due  to  the 

initial  condltlona.  Since  the  remaining  part  of 
the  creep  function  D|(t)  contains  only  retardation 
times  much  larger  than  the  period  of  vibration  we 
may  completely  neglect  the  Influence  of  the  forc¬ 
ing  function.  Eq.  (67)  can  thus  be  simplified  to 

i*  Si  [E  w\  *  £  51  [*  J  «»'(T-*>  ir d#]  • 


8»j|* 

"v 


(74) 


Guided  by  the  solution  of  the  analogous  elas¬ 
tic  problem,  Eqs.  (42)  and  (49a),  lorn's  method  is 
extended  and  a  solution  of  Bq.  (74)  is  sought  la 
the  form 


[5*]  •[«,<’>]  V»,T) 


(75) 


where  *n(B.T)  !•  defined  by  Bq.  (37),  and 


o 


(70) 


[«V<*>]  •  |“(T)  CO*  [*»>(▼)  3  ♦ 

"'ll 

♦  [p(T)/q]  tin  Cq»(T)}/(T,q)  (76«) 

In  Eq.  (76a) 


taken  to  vanish  separately,  flila  results  in 

H')  -  p  *n<T>/q  *  (8U) 

and  In  the  ordinary  differential  equation 

2  ■&  +  ■&-  fo(O)  n  Xs  /  1’  -  0  (81b) 

n 


/(T.q) 


(76b) 


As  In  tha  elastic  problem  all  terms  which  contain 
q  in  the  denominator  are  neglected  in  first 
approximation, 

[wj]  -  ■(▼)  cos  [q»(T)]  vn<M)  (77) 

n 

The  postulated  solution  Eq.  (77)  is  substituted 
into  the  integrodifferential  equation  Eq.  (74). 

For  the  nth  term  this  results  in  tha  aquation 

-  p* X  * (T )  m(T)  cos  [qii>(t)]  V  (I,T)  - 
n  n 

T 

*  q*  j  to(T-»)  **(•)  ■(•)  cos  [qw(s)3  Vn(t,s)ds  • 

o* 


In  evaluating  the  integral  in  Eq.  (78)  we  taka  ad¬ 
vantage  of  the  stipulation  that  g*(T)  covers  the 
part  of  the  relaxation  spectrum  that  consists  of 
larger  relaxation  times,  Iq.  (61b).  ly  integra¬ 
tion  by  parts,  and  by  invoking  the  stipulation  Eq. 
(61b)  the  Integral  is  evaluated  as 

T 

j  g*(T-c)  **(•)  ■(•)  cos  [<r»(s)3  v  (E,s)ds  ~ 
n  n 

o 


■(T)  X*(T)  g,'«» 

q*F) 


sin  [qf«(T)3  Vn(E,T) 


(79) 


Eq.  (81b)  is  satisfied  by 


a(')  -  Eo  exp  [q8g.'(0)T/2p#]/xJ(T)  (82) 

where  E  is  a  constant,  gy  substitution  of  Eqs. 

II 

(81a)  and  (82)  into  Eq.  (77)  we  obtain 

[w$]  -  b(t)  coa  [p  j  X  (s)dsl  V  (X.T)  (83) 

n  o 

It  is  notad  that  in  first  approximation  the 
damping  is  primarily  determined  by  the  larger 
relaxation  times.  The  very  short  relaxation  times 
Influence  the  frequency  of  the  vibration.  For 
fc(0)  ■  0  and  p  ■  q  the  solution  reduces  to  the 
clastic  solution. 

The  complete  solution  W^(E,t)  consists  of  a 

summation  over  n  modes  of  the  type  (83).  The 

constants  X.  are  determined  from  the  initial  con- 
n 

dltloe  on  the  displacement.  For  the  viscoelastic 

problem  tha  displacement  at  t  »  0+  is  obtained 
from  Eq.  (72).  It  is  noticed  that  thla  initial 
value  differs  by  a  multiplicative  constant  N/M 
from  the  dlaplacasMnt  at  t  -  o+  of  the  previously 
considered  elastic  problem.  It  follows  that  the 
viscoelastic  displacement  may  be  written  as 


«i(M).  *(T)  Vn(t,T)cos[p  J  xn(s)ds]  (84) 

q  n  o 


where 

X(T)  -  C(T)  exp  [q8 8b(0)t/p*j  (85) 

and  C(t)  Is  defined  by  Eq.  (56).  Similarly  the 
bond-stress  at  I  ■  1  is  obtained  as 


The  right-hand  aids  of  Eq.  (78)  is  obtained  by 
straightforward  differentiation  of  Eq.  (75).  Aa 
in  solving  the  elastic  problem  Eq.  (78)  is  then 
divided  through  by  Vn(ft,T),  and  the  domain  of 

Interest  is  narrowed  to  the  vicinity  of  X  ■  1  . 

In  view  of  the  limits  (46a)  and  (46b)  the  terms 
containing  dV  (E,T)/dr  and  d  V  (E.'O/dr3  drop 
out.  n  n 

In  the  vicinity  of  E  »  1  Eq.  (78)  reduces  to 

-  p*X*(T)  «(T)  cos  tq«(T)3  * 

-  q  m(t)  X*(t)|9/(0)  sin  [<*i>(t)3  /  4>(t)  - 

•  m(T)  cos  [qu>(T)3-  2  q  m(T)  4i(t)  slnCqw(T)3  - 

-  q  m(t)  tt(T)  sin  [qto(T)3  - 

-  q*  m(T)  d)(T)#  cos  Cq«»(T)3  (80) 

Following  the  general  procedure  of  Itorn's  method^ 
the  coefficients  of  terms  of  ordsr  q  and  q  are 


S)  -  2(Mp*/Nqa)  £  X(T)cos[p  J  X  (s)ds]  (86) 
n  o 

VI.  Discussion  of  the  Eaaulta 

The  axial  shsar  vibrations  of  an  encased 
elastic  cylinder  of  constant  inner  radius  *o  ,  and 

of  elastic  and  viscoelastic  cylinders  of  mono¬ 
tonies  lly  increasing  inner  radii  a(t)  were  studied 
in  this  paper.  The  solutions  of  these  three  prob¬ 
lems  Involved  the  assumption  that  the  uniform 
body-force  distribution  is  suddenly  applied.  The 
Heaviside  unit  function  H(T)  is  used  because  the 
dynamic  effects  are  illustrated  most  effectively 
by  a  suddenly  applied  load.  For  many  practical 
problems  the  load  may,  however,  have  a  finite  rise 
time.  For  such  problems  the  solutions  la  this 
paper  present  upper  limits  on  the  intensity  of  the 
dynamic  effects  that  may  be  expected. 


The  Influence  of  e  finite  riae  cine  can  also 
be  studied  directly  by  modifying  the  present  solu¬ 
tions.  For  the  dynamic  response  of  the  elastic 
cylinder  of  constant  inner  radius  a  this  can  be 

done  very  easily  by  solving  Eq.  (23)  for  arbitrary 
time  dependence  of  An(r).  It  is  easily  shown  tht-t 

the  amplitudes  of  the  vibrations  decrease  with 
increasing  rise  time  of  the  load.  Inclusion  of  a 
finite  rise  time  of  the  load,  for  the  clastic 
cylinder  with  ablating  inner  surface  is  somewhat 
more  troublesome,  since  it  is  then  less  obvious 
what  terms  can  be  ignored  in  Eq.  (49a).  Dynamic 
effects  are,  of  course,  less  pronounced  for  an 
increasing  rise  time.  More  care  must  also  be 
exercised  in  considering  a  gradually  applied  load 
in  the  viscoelastic  cylinder.  If  the  load  is 
suddenly  applied  the  influence  of  the  short  relax¬ 
ation  times  may  be  included  in  the  initial  con¬ 
ditions,  Eq.  (72).  In  this  way  the  damping  influ¬ 
ence  on  the  dynamic  solution  of  the  very  small 
relaxation  times  is  ignored.  If  the  load  is 
gradually  applied  the  Influence  of  the  short 
relaxation  times  further  diminishes.  A  complica¬ 
tion  arises,  however,  when,  depending  on  the  riae 
time,  more  terms  are  needed  In  Eq.  (76a). 

For  the  elastic  grain  of  constant  inner  radius 

R  ■  0  the  shear-bond  overs  tress  at  the  propellant 

casing  Interface,  Eq.  (31),  la  computed  for 

6  ■  0.4  ,  6  ■  0.6  and  6  ■  0.8  .  The  eigenfrequen- 

clas  0  ,  Eq.  (18),  corresponding  to  the  different 
n 

values  of  8,  can  conceivably  ba  obtained  from  Fig. 

1.  The  difference  of  Bessel  functions  In  Eq.  (31) 

la,  however,  very  sensitive  to  small  deviations  in 

0  ,  and  more  accurate  values  of  fi  ^  are  therefore 
n  n 

used  to  compute  the  Bessel  functions  In  Eq.  (31) 

and  in  Eq.  (56). _  It  la  found  that  the  amplitudes 

of  the  modes  of  decrease  rapidly  and  It  la 

sufficient  to  retain  only  the  first  three  modes. 

In  Fig.  2  the  sum  of  the  first  three  modes  and  the 
first  mode  of  the  shear-bond  overstress  are  shown 
for  8  -  0.4  .  The  sum  of  the  first  three  modes  of 
the  shear-bond  overstress  Is  shown  In  Fig.  3  for 
8  »  0.6  and  0  •  0.8  .  It  la  noted  that  for  great¬ 
er  values  of  8  the  amplitudes  are  samller  and  the 
frequencies  are  greater.  The  relatively  soft 
material  of  the  cylinder  Is  defined  by  q  -  5.103, 
Eq.  (9c). 

The  shear-bond  overstress  at  R  »  1  for  the 
ablating  elastic  grain  Is  examined  for  various 
ablation  rates.  The  Influence  of  the  ablation 
function  o(t)  Is  shown  by  considering  the  func¬ 
tions 

Qf(r)  -  1  +Xt  ,  where  X  -  (1  /S)  -  1  (87a) 

and 

er(T)  •  (1-xt)"^  ,  where  h  ■  1-B8  (87b) 

For  8  ■  0.4  the  ablation  functions  (87a)  and 
(87b)  are  shown  In  Fig.  4. 

From  the  expression  tor  1^(8, t),  Eq.  (57),  It 
Is  noted  that  the  amplitudes  of  the  modes 
decrease  and  the  frequencies  Increase  towards 
burnout.  Me  define  the  frequency  function 
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where  the  constan 

ts  c, 

.2,3 

and  d,^^ 

are 

shown  In 

Eqs .  (41a, b,c).  By  subsLitu'ion  of  Eqs.  (87a), 
(87b)  and  (89)  into  Eq.  (88)  the  frequency  func¬ 
tions  are  evaluated  as 


*  .  -  d  t  -  [(c  -d  )/(1-«)]  t n  (1-T)  (90a) 

m  n  n  n 

♦  ,  •-  [e  /(I -r* )Hi -KT)%  +  6]*  ♦  c  (l+8)/(i-8)  ♦ 

n<  n  n 

+  2d  e[(1-KT)*  -  1  3/(1  -ft*)  -  (90b) 

n 

-  t28a(cii-dti)/(1-e*)]  la  {t(1  -Kt)4 -  B]/(1  .B>} 


The  equations  (90a)  and  (90b)  Indicate  that  the 
frequencies  Increase  rapidly  towards  burnout  as  t 
approaches  unity.  From  Eq.  (56)  It  Is  noted,  how¬ 
ever,  that  the  amplitudes  approach  sero  as  T 
approaches  unity. 

The  frequency  functions  j  for  the  first  mode 

arc  shown  In  Fig.  5.  The  frequency  functions  are 
determined  up  to  t  ■  0.95.  It  Is  noticed  thst 
linear  ablation  according  to  Eq.  (87a)  results  In 
a  much  more  rapid  Increase  of  the  frequencies. 

The  frequency  functions  ♦  .  for  higher  modes  arc 
shown  In  Fig.  6.  n 

The  observations  on  asq>lltudes  and  frequencies 
are  further  illustrated  by  Figures  7,  8,  9  end  10. 
The  figures  7  end  9  show  the  first  mode  of  Eq. 
(57),  respectively  near  t  •  0.5  and  t  •  0.9  .  The 
sum  of  the  first  three  modes  of  the  shear-bond 
stress,  Eq.  (57),  is  shown  for  t  -  0.5  and  t  ■  0.9 
In  Fig.  8  and  Fig.  10  respectively. 

The  results  presented  In  this  paper  show  that 
the  frequencies  of  axial  shear  vibrations  Increase 
significantly  during  ablation.  This  effect  may 
Influence  the  structural  Integrity  of  the 
propellant-casing  system. 
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Fig.  I  Frequency  coefficients  fi n  versus  $. 
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Fig.  2  Shear- bond  overstress ;  no  ablation  fi  s  0.4 


Fig.  3  Shear -bond  overstress;  no  ablation  £*o.e  and 


Fig. 4  Ablation  function  o(t)  versus  t  for  P  -  0.4 


Fig.  6  Frequency  functions  ^ni  versus  t  for  higher  modes. 
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Fig. 7  The  first  mode  of  the  shear -bond  overstress 
for  various  ablation  functions  ;  x  *  0.5  x'. 
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Fig. 8  Sheor-bond  overstress  for  T«o.5■♦■T, 
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Fig.  9  First  mode  of  the  sheor-bond  overstress  for 


Fig.  10  Sheor-bond  overstress  for  t  *  o.t  ♦  xf 


